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Abstract. Quantum vacuum energy (Casimir energy) is reviewed for a mathematical 



I audience as a topic in spectral theory. Then some one-dimensional systems are solved exactly. 



in terms of closed classical paths and periodic orbits. The relations among local spectral 
densities, energy densities, global eigenvalue densities, and total energies are demonstrated. 
I This material provides background and motivation for the treatment of higher-dimensional 

systems (self-adjoint second-order partial differential operators) by semiclassical approxi- 
mation and other methods. 
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1 Introduction 

I Vacuum energy is a concept arising in quantum field theory, with observable consequences of 

J> ' considerable current interest in physics [13l [52l [501 [71 HI] . Here, however, I treat it as a purely 

mathematical problem, an underdeveloped aspect of the spectral theory of self-adjoint second- 
(yQ ! order differential operators. 

CN ' Vacuum energy is related to the oscillatory terms in the spectral density, which are associated 

' with the periodic orbits in the classical mechanics determined by the operator as Hamiltonian 

, [H [21 [T5l [36l [T9l [T6l [9] . The information concerning these oscillations is lost from the much- 

Q I studied short-time expansion of the heat kernel [SH [Ml [351 1321 M US] > but some of it remains 

in the asymptotics of another integral kernel, the cylinder kernel, which includes the vacuum 
energy directly [HI [5l [371 IMl [Ml [27] . The local density of vacuum energy is related analogously 
to a local spectral density, associated with classical paths that are closed but not necessarily 
' periodic. Thus the vacuum energy and its density are probes of the detailed geometry associated 

with the operator. 

The connection between closed classical paths and Casimir energy is implicit in the two 
papers of Balian and Duplantier [3], which treat those two aspects of the electromagnetic field 
separately. In recent years it has become a serious tool for the calculation of vacuum energies 
[351 [571 B9l [271 [5^ [55l [56l [39] . In one spatial dimension it reduces to the classical method of 
images, which has been applied to Casimir calculations since at least 1969 |11] . 

The main purpose of this article is to demonstrate in detail, for the simplest one-dimensional 
models, the connections among closed orbits, spectral densities, and vacuum energy densities, 
as well as the corresponding global quantities. Sections [3] and [H are based on research notes 
that I have used for some years with my students and collaborators; the fact that they have not 
been published and available for citation has become increasingly inconvenient. First, however, 
in Section [2] I present the basics of vacuum energy in a more general context to a mathematical 



*This paper is a contribution to the Proceedings of the 2007 Midwest Geometry Conference in honor of 
Thomas P. Branson. The full collection is available at |http: / /www. emis.de/journals / SIGMA/MGC2007.html , 



2 



S.A. Fulling 



audience. This material is based on talks given at the Workshop on Semiclassical Approximation 
and Vacuum Energy at Texas A&M University in January, 2005, and the Workshop on Spectral 
Theory and Its Applications at the Isaac Newton Institute, University of Cambridge, in July, 
2006. I expect to publish elsewhere [30] a more complete discussion of the physics of the subject, 
still for a mathematical audience, along with a review of the recent work on vacuum energy in 
quantum graphs |59j . 



2 Vacuum energy (and energy density) in general 
2.1 Spectral theory 

Let if be a second-order, elliptic, self-adjoint partial differential operator, on scalar functions, 
in a d-dimensional region Q. The prototype situation is a billiard: 

H = -v^ n c R'^, 



with boundary conditions that make the operator self-adjoint (such as the Dirichlet condition, 
u = on di}). The treatment can be generalized to the electromagnetic field (vector functions) - 
which is the case of greatest physical interest - and to other boundary conditions, to Riemannian 
manifolds (Laplace-Beltrami operators), potentials {H = — V'^ + y(x)), applied magnetic fields 
(gauge-invariant Laplacians on vector bundles), etc. For simplicity, assume that the spectrum 
of H is nonnegative and that if is an eigenvalue, then the eigenfunction is constant. Of course, 
for precise theorems some smoothness hypotheses on the boundary, potential, etc. are needed. 

For the moment let us assume that the spectrum is discrete, as will be so if is a com- 
pact billiard. In this case a finite total vacuum energy is expected. We review the spectral 
decomposition and functional calculus. Let (pn be the normalized eigenfunctions: 



Hipn = \n'^n, lbn|P= / \pnix)\'^ dx = 1. 

Jn 

Define u;„ = ^/X^. Functions of the operator H are defined by 

CO „ 

f{H)u = 'S2f{Xn){'Pn,u)ipn, {ipn,u) = / (pn{x)u{x) dx . 

n=l 

At least formally, f{H) is given by an integral kernel: 

„ oo 

f{H)u{x)= I G{x,y)u{y)dy, G{x,y) = f{\n)pn{x)ipn{y). 

n=l 

In general G is a distribution, but if / is sufficiently rapidly decreasing, G is a smooth function, 
and the trace is defined: 



„ CO 

TrG= / G(x,x)dx = V/(A„) = Tr/(ii). 

n=l 



The prototype is the heat kernel, G{x,y) = K{t,x,y), corresponding to the parametrized 
function /t(A) = e"*"^. Then u{t,x) = ft{H)uo solves 

du 

— = -Hu, u{0,x) = uo{x). 
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It is well known [351 [El [321 [IS] that K for a billiard has the asymptotic expansion 

oo oo 

TrK = J2 e"*^" ~ bst^"^^'^'^ (2-1) 

n=l s=0 

where each term has a global geometrical significance - for instance, 6o is proportional to 
the volume of The inverse Laplace transform of the leading term in Th: K gives the leading 
behavior at large A of the density of eigenvalues, and the higher-order terms correspond similarly 
to lower-order corrections to the eigenvalue distribution, on the average [T^l HHl [53] . (When V{x) 
is a confining potential [361 [H [9], i/ = — + V may have discrete spectrum even though its 
spatial domain, 0, is not compact. In such a case the form of the asymptotic expansion of Tr/C, 
and the resulting asymptotic expansion of the eigenvalue distribution, may be rather different. 
For example, the eigenvalues of a one-dimensional harmonic oscillator are evenly spaced in A, 
whereas those of a one-dimensional billiard are evenly spaced in w.) 

The cylinder kernel (also called Poisson kernel), T{t,x,y), corresponds to the parametrized 
function /t(A) = e~*^. That is, ft{H)uQ is the solution of 

= Hu, u(0, x) = uoix), 

that is well-behaved as i — > +oo. (T is a boundary value of a derivative of the fundamental 
solution of the elliptic operator H — ^ in 17 x R. It is the ("imaginary-time") analytic conti- 
nuation of the time derivative of the Wightman function, a certain fundamental solution of the 
hyperbolic operator H + ^ .) We have 

oo „ oo 

T(t, x,y) = Y, e~''^-^n{x)^n{y): Tt T = T{t, X, x)dx = Y, 

n=l •'^ n=l 



As t I one has the asymptotics [Ml [23 [1] 

oo oo 

TrT- J^e,t-'^+^+ /.i-'^+Mnt, (2.2) 

s=0 s=d+l 
s~d odd 

as described in the following theorem. (For a more complete, but succinct, statement of the 
connection with Riesz means, see [27]. See also [HI [201 [211 [2S]-) 

Theorem 1. The traces of the heat kernel and the cylinder kernel of a positive, self- adjoint, 
second-order linear differential operator in dimension d have the asymptotic expansions (|2.1|) 
and (j2.2|) . and precisely parallel expansions hold for the local (untraced) diagonal values of those 
kernels. The hg are proportional to coefficients in the high-frequency asymptotics of Riesz means 
of N (or P) with respect to A. The Cg and fs are proportional to coefficients in the asymptotics 
of Riesz means with respect to uj. If d — s is even or positive, 

es=7r-'^^2'^~'T{{d-s + l)/2)bs. 

If d — s is odd and negative, 

(^_-Y'^(s-d+l)/22d-s+l 

" ^T{{s-d+l)/2) 
hut Cs is undetermined by the br- 
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2.2 Vacuum energy 

We can now define the vacuum energy as the coefficient with d — s = —1, 

E = -\ei+d- 
Formally, E is the "finite part" of 



(2.3) 



1 1 d X - 



n=l 



t=0 



(When fi^d / 0, -E is actually defined only modulo a multiple of /i+rf , because of the scale 
ambiguity in the argument of the logarithm. That complication will not arise in the problems 
studied in this paper.) 

The prototype example is 



d nvr /mrx 

n={0,L), if=-_; uJn = —, V3„,(x) = Sm ^ 

One can evaluate the cylinder kernel directly from the spectral decomposition as 



(2.4) 



T{t,x,y) = j^sin 



n=l 



mrx 



sm 



or by the method of images as a sum over classical paths, 



T{t,x,y) 



t 



E 



1 



Af=-oo 

Either sum can be evaluated in closed form as 



(x - y - 2iVL)2 + t2 (a; + y _ 2ArL)2 + t2 



(2.5) 



T{t,x,y) 



It follows that 



1 

2L 



sinh(7rt/L) 



sinh(7rt/L) 



cosh(7rt/L) — cos(7r(x — y)/L) cosh(7ri/L) — cos(7r(x + y)/L) 



(2.6) 



^ „ 1 sinh(7rt/L) 1 L 1 irt ^, 

TrT = , , , \ + + 0(t^) 

2 cosh(7rt/L) - 1 2 vri 2 12L ^ ^ 

Thus the energy, the 0{t) term times — ^, is 

TT 

E -- 



24L 



(2.7) 



This formula has been known for many years (e.g., [5]). 

Another simple example, or class of them, involves a vector bundle over the circle, coordina- 
tized as = (0, L) [ISIEH]- Again H = -d^/dx"^, but now the field is "twisted" so that 



u{L) = e^%(0), u'{L) = e*%'(0). 
The eigenfunctions and eigenvalues are 

J(2-Kn+B)x/L) t„ ^ ry\ , ,± _ ^TTJ ± 



iPn{x) = e* 



(n G Z), 



L 



(j G N or 



(2.8) 



(2.9) 
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One can then show (see Section [3]) that 



Ea 



IT 



TT 

UL 



6- + 3 

TT \TT 



(2.10) 



(i?2 is a BernoulH polynomial.) What is intriguing about this example is its dependence on the 
parameter 9; Eg can be positive as well as negative. Its extreme values are Eq = — and 
E'jT = + ; these are the cases where the two eigenvalue sequences coincide, so the gaps in the 
spectrum are largest. Eg passes through for ^ 0.42, which is close to the point tt/2 where 
the eigenvalues are equally spaced. We observe that for an individual eigenvalue sequence, the 
sign and magnitude of its contribution to the energy are determined by the phase of the spectral 
oscillation. Here that phase is controlled hy 6; in higher-dimensional systems the phase of the 
oscillations associated with a periodic classical orbit is controlled by the famous [Kramers- 
Morse-Keller-Gutzwiller-jMaslov index (see, e.g., [E]), whose significance for vacuum energy 
has only begun to be explored [551 l56l [23]. 

2.3 Vacuum energy density 

For present purposes the energy density can be defined simply by leaving out the integration 
over X in the cylinder trace and proceeding as before. Let P{\, x, y) be the integral kernel of 
the orthogonal projection onto the part of L^(r2) corresponding to spectrum less than or equal 
to A; since 



/•oo 

K{t,x,y) = / e-'^dP{X,x,y), 
Jo 



we can also define P as the inverse Laplace transform of the heat kernel (the exact heat kernel, 
not its small-t series). Then 



T{t,x,y) = I e 

'0 



~'^dP{X,x,y), 



and, as mentioned in Theorem [H its trace has an expansion completely analogous to (12. 2p with 
x-dependent coefficients: 



r(t,x,x)~5^e,(x)t-'^+^+ J2 fsix)t-'''-nnt. 



s=0 



s=d+l 
s—d odd 



We can now define the energy density as 
E{x) = -iei+d(x). 

In the quantum field theory, E(x) (known in that context as Too(x) defined with = :j) is 
formally the finite part of 



-1 +uHn 



When the spectrum is discrete, 

P{\x,y) = ^ Lpn{x)<fn{v), 



A„<A 
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and the integral of P{X,x,x) over equals N{X), the number of eigenvalues less than or equal 
to A. But one reason for the importance of the energy density is that it remains meaningful 
when H has some continuous spectrum. If the spectrum is absolutely continuous, 

a{uj,x) = —P{u>'^,x,x) 
dio 

exists (as a Radon-Nikodym derivative) and constitutes a local spectral density. (It is a density 
in two senses - with respect to both lo and x.) 

Here, also, there is a prototypical example, the half-line with a Dirichlet endpoint: 

J7 = (0,oo), H = --^, n(0) = 0. (2.11) 
The eigenfunction expansion is the Fourier sine transform, equivalent to the projector 



f 2 

P{X,x,y) = / — sm{kx) sm{ky) dk. 
Jo 



/o 

The solution for the cylinder kernel by the method of images is now a simpler analogue of (j2.5p : 

1 1 



T{t,x,y) = - 



{x-yf + t"^ (x + y)2 + t2 
Thus 

2k 



(2.12) 



1 f / f \' 



so 



Recall that in any one-dimensional billiard the untraced diagonal value of the heat kernel is 
simply 

K{t,x,x) ~ [ATTtr^/^ + 0{t°°) 

(for fixed x in the interior of $7) regardless of boundary conditions. In contrast, we see from (j2.6p . 
(|2.10p . and (j2.13p that T{t,x,x) is sensitive to the global geometry - the length of the inter- 
val, the nature of the boundary conditions, the structure (angle 6) of the vector bundle. (In 
case ()2.10p . E{x) is independent of x and equals E/L.) Therefore, T{t,x,x), E{x), and E are 
interesting objects to study from the point of view of inverse problems and other aspects of 
spectral geometry. 

For the finite interval (12.41). one finds 



This function is graphed in Fig. [TJ Now 

TT 2 /"^^^ 1 



with a similar expansion as x — > L. Therefore, E{x) displays the boundary behavior (|2.14p 
in addition to a spatially homogeneous term. But now notice something strange: The "total" 
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Figure 1. Boundary energy density for il — (0, 1). Figure 2. Regularized energy density E{t, x) for 

n = (0,oo) and t = 0.001. 



energy E, (j2.7p . is equal to the integral over Q of the spatially homogeneous term only, while 
the integral of the boundary terms diverges! In physical terms, the renormalized energy is not 
equal to the integral of the renormalized energy density. The physical significance of this fact is 
beyond the scope of this article. (In fact, it is still controversial.) Mathematically, however, it 
is just an instance of nonuniform convergence: The limit t — > cannot be interchanged with the 
limit X — > 0, and hence with the integration over x. 

To examine this phenomenon more closely, it is convenient to return to the half-line (j2.1ip . 
Since the latter is spatially infinite, one would not necessarily expect its total energy to be 
defined; however, we notice that the integral of (j2.14p converges at infinity. Let us keep t 
positive (but small) and consider the regularized vacuum energy density 

E{t, x) = -- -rnTit, x,x)-—^ 



2dt ' ' ' ' 7rt2 27r (t2 + 42;2)2 " 

(See Fig. [2l noting the compressed vertical scale. The subtracted term is the contribution of the 
first term in p.l2p (or the first term in ()2.13p ). identified with the ubiquitous but unobservable 
vacuum energy of infinite empty space. Only the second term in (j2.12p is of interest.) The limit 
t — > reproduces the renormalized energy density (I2.14p . which would imply an infinite total 
energy if integrated. However, 



/•oo 

E{t) = / E{t, x)dx = for ah t > 0, 
Jo 



which implies a vanishing total energy if one now takes t — > 0. This disappearing divergence 
is a general property of energy densities that behave like x~'^ near a boundary. (It has been 
confirmed for the corners of a rectangle or rectangular parallelepiped [23].) By dimensional 
analysis any correlate of such a term in E{t) must be proportional to t~^ and hence must come 
from a term /olnt in (12. 2p . but no such term can exist. (It would necessarily match a Int term 
in the heat kernel [14 1 1331 [6ll26|. which does not exist for the class of operators considered here.) 
For the finite interval, it will be shown in Section S] that the regularized total energy does not 
vanish, but does approach the renormalized energy ()2.7p as t — > 0. 



2.4 But what about the zeta function? 

Many mathematically inclined researchers on vacuum energy have defined it in terms of zeta 
functions. For example, (|2.7p is often expressed as 

2 ^ L 2L^^ ' 24L' 

n=l 
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where C is the original Riemann zeta function. Therefore, a comment is needed upon the relation 
of generalized zeta functions to the approach expounded here. 

Let s be a complex parameter, and consider the operator function fsiH) = H^^ . Then the 
zeta function for H is defined by C{s,H) = Tr fs{H) and extended by analytic continuation to 
values of s for which 

oo 
n=l 

does not converge. Note that 

C{s,H) = C{2s,Vh). (2.16) 
On the other hand, these zeta functions are related to our integral kernels by 

/•oo 

Jo 

and a similar equation with K{t, H). It follows from (j2.ip and (j2.2p that 6„ and are residues 
at poles ofT{s)C,{s,H) (at s = ^{d — n)) and r{s)({s,y/H) (at s = d — n), respectively. So when 
there are no logarithmic terms, (I2.16P implies 

r(^)"\„ = ^r(d-nrV 

(This observation is analogous to one by Gilkey [31] concerning higher-order operators.) Now 
when d—n is odd and negative, T{d—n) has a pole where T(^^{d—n)) does not; the information in 
the corresponding e„ is thereby expunged from the heat-kernel expansion, and that is how (j2.ip 
contains less information than (j2.2p . In that case, e„ ceases to be a residue of the zeta function 
and becomes a value of zeta at a regular point - a more subtle object to calculate. (Logarithmic 
terms give rise to coinciding poles of ^ and F.) 



3 Vacuum energy in a vector bundle from periodic orbits 

Let us concentrate now on the model studied in [28] and summarized above. H = —-^ acts in 
L^(0, L) with the modified periodic boundary conditions (j2.8p . representing a nontrivial holono- 
my in a line bundle over the circle. For definiteness take 6 E [0, 27r). (A gauge transformation, 
(p{x) = e^^^^^ip{x), converts this problem to "Bohm-Aharonov" form with ordinary periodic 
boundary conditions and a nonzero, but pure gauge, vector potential.) Eigenfunctions must be 
proportional to e^^^^ with kj = (27rj + d)/ L, j ^ Z; the eigenvalues for positive and negative k 
then form the two sequences ()2.9p . 

Because of the spatial homogeneity of the model, there is no distinction, except a factor L, 
between total energy and local energy density. Although vacuum energy is barely mentioned 
in [28], (j2.10p is equivalent to either of the two main results of that paper, the cylinder kernel for 
H and its small-t diagonal expansion [28^ (9) and (13)], and the first nonlocal Riesz mean of the 
eigenvalue density with respect to frequency iv [281 equation above (30)]. Here I rederive (I2.10p 
in a third way, in the framework of periodic-orbit theory [H \T5\ [T9] (but using the exact formula 
for the Green function, not a stationary-phase approximation). 

The Green function (resolvent kernel) for the time-independent Schrodinger equation in 
infinite one-dimensional space is 

Goo(a;^x,y) = ;^e*-l^-^l. 
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(By convention, for llP' on the positive real axis is defined as the limit from above 
(J + ie), or, equivalently, by the outgoing radiation condition. Also, let us adopt the convention 
{Hx — \)G{X,x,y) = +6{x — y), where much of the physical literature has the opposite sign.) 
The Green function for the bundle is constructed as a sum over all classical paths connecting y 
to X in the covering space; in this simple one-dimensional problem this construction reduces to 
the traditional "method of images". When y = x the paths become periodic orbits, and there 
is one for each n G Z, with length \n\L. Since there are no boundaries to reflect the paths, in 
this model there is no distinction between periodic orbits and more general closed paths (with 
different initial and final velocity), in keeping with the previous observation that energy and 
energy density are the same thing. The Green function G = G^^q is 

oo 

G{uj^,x,y)= Goo{oj\x,y + nL)e'^'. 

n=— oo 

(The sum converges distributionally and has the desired twisted-periodicity property.) 

The starting point of periodic-orbit theory is that the density of eigenvalues, as a function 
of A, is 



= -lmTrG(A) 



The trace is an integration over x, amounting here to a factor L. It is more convenient to 
work either with the density with respect to w = \/A, which carries an additional factor 2u;, or 
with the eigenvalue counting function N , which is the same quantity whether A or oo is used as 
independent variable. We have 

oo 

ImG(u;^x,x) =Im V J-eH^l^e^"^ 

n=— oo 
1 1 °° 

= \ y [cos(ti;nL + nO) + cos{ujnL — nO)]. 

n=l 

(One can identify tij|?i|L as the action of the nth periodic orbit. There is no Maslov index in 
this problem, but there is a phase shift nO from the nontrivial holonomy.) So the integrated 
eigenvalue density is (for > 0) 

L f Llo 1 1 
N(uj) = — / 2uJ duJ ImGiuj'^ , x,x) = 1 —[sm(u!nL + nO) + siniujnL — nO)] 

•' n=l 
= N^,{uj)+Nosc{i0). (3.1) 

Here Ns;v{^)/L = w/vr (which came from the orbit of zero length) is the density of states per 
unit length that would exist in infinite space, and Nosc{^) describes the bunching of spectrum 
caused by the existence of closed orbits. (The contribution of orbit n to (j3.ip includes an to- 
independent constant of integration, — sin(n0)/n7r, which cancels with the corresponding term 
from orbit — n.) 

Before turning to the vacuum energy, let's digress to see how (jS.ip reproduces the known 
eigenvalues. It is known [Ml 1.441.1, 9.627.1] that for < z < 27r 

oo 

Vsmnz z 
= TT = — vrisi 
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where Bi is the first Bernoulli polynomial; elsewhere, the sum defines a 27r-periodic sawtooth 
function, vanishing at the discontinuity points z = 2j7r. Thus, for a; > 0, 

Noscico) = ;^[/(vr - - 0) + /{tt - loL + 9)] 
Zn 

where /(C) is the 27r-periodic extension of the function equal to C when — vr < ( < ir and equal 
to at the endpoints. Hence /(C) jumps by — 27r at each odd multiple of vr. In more transparent 
terms, the eigenvalue density is (since A'av and A'osc vanish for u; < 0) 

j=l ^ ^ j=0 ^ ^ 

in precise agreement with the eigenvalues calculated already by elementary means. Of course, 
in higher-dimensional problems such exact results are not to be expected. 

The renormalized vacuum energy is the contribution of A'^sc ^o ^ ujp{uj) duj. Returning 
to (j3.ip . and considering the local energy density, we have (ignoring analytical technicalities for 
the moment) 

E{x) = I — uj N^g^{u>) du! = — > / [cos{u!nL + n9) + cos{ijL!nL — n6)]uj du). 
Jo 2L 2tt ^ Jq 

Consider just one term, and integrate to a finite upper limit: 

/ cosiujnL + nO) ujduj = — - smi^nL + nO) + - — — -r cos{VtnL + n9) — - — — -75- cos(n0). 
Jq nL [nL)''' [nLy 

As — > 00, the first two terms oscillate with zero mean. If we assume for the moment that 
they can be ignored, we have 

^ cos(n6l) (3.2) 



27r(nL)2 



as the contribution to the energy density from one of the two periodic orbits of length nL; that 
from the other orbit comes out the same. 

Again from [Ml 1.443.3, 9.627.2], for < 6" < 27r 



00 



^ cosn^ 7r2 ^9 9"^ 2„ M A 

71=1 ^ ^ 

where B2 is the second Bernoulli polynomial. The total energy density is thus 

^ ^ R M 

^^-) = -T^^^ (2^J=T' 

Eq given by (I2.10|) . as was to be verified. 

It is noteworthy that the vacuum energy (j3.2p associated with a single spectral oscillation 
depends critically, in sign as well as magnitude, on the phase n9 of the oscillation. Algebraically, 
(j3.2p comes entirely from the lower limit of the integration over and an effective lower cutoff 
on that integration would appear to change the phase. This is a matter of great concern for the 
extendability of the theory, since in more general circumstances the Gutzwiller spectral oscilla- 
tions arise from stationary-phase approximations that are not justified at low frequency. It has 
often been observed, however, that periodic-orbit calculations reproduce low-lying eigenvalues 
more accurately than they have any right to do. 
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The problem at high frequency, in contrast, is more apparent than real. The integral over to 
can be defined by Riesz-Cesaro summation of order 2: A lengthy exercise in integration by parts 
shows that 



/■^ / a;\2 1 

/ ( 1 ~ o ) cos(a;n-L + nO) oj duj = — - — —-r cos{n6) + 0(0 



Taking the limit Q ^ oo now yields (j3.2|) . Clearly, the interchange of integration and summation 
in the calculation of E{x) is now also justifiable. 

Another approach leading to the same conclusion is Abel summation, which amounts to 
finding the contribution of each spectral oscillation to the cylinder kernel. We are after the 
small-t behavior of 



/■oo ^ ^ roo 

/ -ujoiuj^e^"^^ duj = -/ piuj^e^"^^ duj. 

Jo 2 ^ '^dtJo ^ 

Using the integral 







cosfaw — b)e '^^ duj = -75 ^ cos 6+ sin 6, (3.3) 

+ + 



one finds that the contribution of the two orbits with length 7iL to p{^)g du is 
2L t 



vr t^ + {nLy 



cos{n6). 



(The two contributions are not equal; their sine terms cancel. But because the sine term in (13. 3p 
has no term of order exactly in its small-t expansion, it would make no contribution to (j3.2p 
even before the pairwise cancellation.) The sum over n is now absolutely convergent, and at 
t = the expected formulas (j3.2|) and ()2.10p emerge. Alternatively, the series at finite t can be 
summed in closed form |341 1.445.2], and with the help of some identities for hyperbolic functions 
the result is shown equal to the diagonal value of the cylinder kernel as found previously [281 
(12)]. 

If 9 is written \9\, ()2.10p is valid over the interval — 27r < 9 < 27r. Overall, Eg is a continuous 
27r-periodic function with cusps at the integer multiples of 2ir; it is symmetric under reflection 
about any integer multiple of vr. We have already commented in Section [2] about the significance 
of its extreme values and its zero. 



4 Boundary vacuum energy from closed and periodic orbits 

We consider a finite interval with either a Dirichlet or a Neumann boundary condition at each 
end. Thus H = — ^ acts in L^{0, L) on the domain defined by 

(0) = 0, 4^-''^ (L) = 0, where /, r G {0, 1}. 
For / = 1, the (unnormalized) eigenfunctions are sm.{LOjx) with 

ujj = j-, je Z+, if r = 1; cvj = ^ (j + ^) > J G N, if r = 0. 
For / = 0, they are cos{ujjx) with j E N and 
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The Green function can be constructed from Goo by the method of images, or, equivalently, 
as a sum over reflected paths as in |57[ Fig. 1]. (I suppress the arguments and x after their 
first appearance.) 

G{oj\x, y) = Goo(y) + (-l)'Goo(-y) + (-irGoo(2L - y) + (-l)'+'^Goo(2L + y) 
+ (-l)'+'-Goo(-2L + y) + (-l)2'+-Goo(-2L - y) 
+ (-l)'+2'-G„,(4L -y) + (-l)2'+2r^^(4^ + y) + • • • 



£(_l)"a+0G'^(_2nL + y)+ ^(-l)'+"('+'')Goo(-2nL - y) 

n=0 n=0 

oo oo 

+ 2](-l)-'+"^'+'')Goo(2nL - y) + ^(-l)"('+^)Goo(2nL + y) 



71=1 71 = 1 

OO . OO 



l\n{l+r) M2nL+x-y\ , _J_Sr^ (-^\l+n[l+r) iu{2nL+x+y) 



2u;^' ' 2u; 

n=0 71=0 

oo . oo 

_|_ _L ''^(^_Yyl+n(l+r) ^iw{2nL-x-y) _j_ J_ ^|-_-|^-jn(«+r)gja;(27iL-a;+3/) 

n=l n=l 

(cf. |571 (2.16)]). When y = x, the first and fourth sums correspond to periodic paths (including 
the direct path of zero length), but the second and third sums come from paths that are closed 
but not periodic (i.e., the "particle" has hit the boundary an odd number of times and returned 
to its starting point with reversed velocity). The n = term in the first sum is the direct path; 
n = in the second term and n = 1 in the third (which will become n = —1 in the next step) 
are the "short" paths that bounce off the boundary only once. 
We have (always with the understanding that uj > 0) 

vrcr = 2u!lmG{uj'^ ,x,x) 

oo oo 

= X^C-l)"^'^"^ cos(2wnL) + ^(-l)'+«(^+0 cos(2u;(nL + x)) 

n=0 n=0 

oo oo 

+ ^(-l)-^+'^('+^) cos(2cj(nL - x)) + ^(-1)'^('+^) cos(2cjnL) 

n=l 71=1 
oo oo 

= 1 + 2 cos(2u;nL) + ^ (_i)«+n(«+r) cos(2w(x + nL)) 

71=1 71.= — OO 

= 7r((Tav + fpcr + <7bdry) = 7r((7av + CFosc)- 

Because the series are not absolutely convergent, we must be wary of cancelling contributions 
from different values of n. For the most part, we want to consider to the latest possible moment 
the contribution from each orbit separately, and the summation is just a reminder of which 
values of n occur. We do pair n with —n or — (n + 1) when their contributions are manifestly 
identical or parallel. We have just done so for the first and fourth terms in a, which became 
o'av + o'per- In the second and third terms, Ubdry; it will be natural to pair n ( > 0) with — (n + 1), 
since that combines paths that have the same number of reflections, 2n + 1 (again consult [571 
Fig. 1], or just look at the exponents of (—1) in the formula). At fixed x, orbits n and — (n + 1) 
do not have the same length, but when integrated over x those two sets of orbits collectively 
occupy the same part of the length spectrum. 

Let us compare with the case of a single boundary (Hilbert space L'^{0, oo), boundary condi- 
tion «(^-')(0) =0). In that case only the direct path and the short path off the boundary exist. 

Thus cTper = 0, (Jav is Unchanged, and o"bdry = ^ cos(2a;x) has only one term. 
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There are now two routes to follow: We can work with the local spectral density a (integrating 
with respect to to, summing over n, calculating the local energy density, etc.), or we can integrate 
over X first to get the global eigenvalue density. (Unfortunately, "density" is ambiguous in this 
context, as previously remarked.) 

Let's examine the global situation first. The eigenvalue density is 



p{uj) = / a{u;, x) dx = pweyl + Pper + Pbdry , 

Jo 



where 



PWeyi = / <ya.wdx = -, Pp^r = — V'(-l)"('+'') cos(2a;nL), 

n=l 



Pbdry 



— 2^ [sin(2a;L(n + 1)) - sin(2u;Ln)]. 



27r ^-^ UJ 

n=— oo 



In keeping with previous remarks, we refuse for the moment to take advantage of the possibility 
of "telescoping" pbdry when / + r is even, but we can combine the positive and negative parts of 

the series now. In the summand — '— [sin(26<jL(n + 1)) — sin(2ti;Ln)], replace n by — (n + 1): 

(-l\-in+m+r) 

^ ' — [sin(2a;L(-n)) - sin(-2a;L(n + 1))] 



UJ 



(_1\n{l+r) 

(^_Y^l+r\ 1 [sin(2wL(n + 1)) - sin(2a;Ln)], 

UJ 



which is the original summand except for the initial sign. Therefore, if Z + r is odd (the two 
boundary conditions are different), pbdry = 0; but if / + r is even. 



(-1)' °° 1 

Pbdry = / — [sin(2a;L(n + 1)) — sm.{2ujLn] 



vr — ' UJ 

n.=0 

(The local counterpart of this observation is that the simultaneous transformation — (n+1), 
X I— > -L — X leaves the summand invariant up to the sign. Together with the sign (-1)', this 
indicates that the effects that are localized near the boundaries are equal and opposite for 
Dirichlet and Neumann boundaries.) 

The corresponding equation for the half-line with a single boundary is formally 

Pbdry = ^ / cos(2u;x) dx = ^-^5{uj). (4.1) 

Jo 4 

Of course pav is infinite in that case because of the infinite volume and continuous spectrum, 
but Pbdry ! being associated with the boundary, is a localized, finite contribution that can be 
studied separately. 

Remark. The distributional integral (14. ip has been thoroughly studied in [22]. The second 
equality in (|4.ip is correct under the convention that 

/■oo poo 

f{0)= / 6iuj)fiu;)duj= / 6{uj)H{uj)fiuj)du;, (4.2) 



where H is the unit step function. An alternative convention is that a delta function at an 
endpoint of an interval of integration yields only half the value of the test function at that 
point; otherwise put, the integration is extended over all a; E R and the integrand is interpreted 
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according to the rule 6{u>)0{u>) = ^6{oj). The choice is somewhat arbitrary [22], but (j4.2p has the 
advantage that an eigenvalue ujq always corresponds in p to a unit-normalized delta function, 
5{uj — loq), even when = 0. 

In a more general problem, the integrals over lo would need to be evaluated in a stationary- 
phase approximation. It would then be argued that only periodic orbits contribute, so /Obdry 
would be set to zero. A formal justification for this approximation is that the neglected terms 
are of higher order in Planck's constant when nonrelativistic quantum-mechanical units are 
used in H. In relativistic vacuum-energy calculations in natural units, the manifestation of 
this observation is that Pbdry is suppressed relative to pav + Ppcr by a factor 1/loL, which is 
small in the high-frequency regime where WKB-type asymptotics would be valid. A less drastic 
approximation is to keep only the two short orbits: 

(-1)^ sin(2a;L) 

/5bdry ~ • 

TT W 

This is plausible because the boundary effect should come only from points close to the boundary 
(cf. the local calculations below and the single-boundary equation above). 

The eigenvalue counting function N{uj) is zero for uj < and p for w > 0. Therefore (for 
w > 0), 

^Weyi = — , ^pcr = - sin(2a;nL), 

n=l 

' (-ly ^ r sin(2a;L(n + 1)) - sin(2a;Ln) , 
> / du! if / + r IS even. 



if / + r is odd. 



We contemplate each of these in turn. 
A^Weyi is exactly as expected. 
Ul + ris odd, by [Ml 1.441.3] we have 

A^pcr = - V sin(2a;nL) = - — /(2wL) 

TT n ZTT 

n=l 

with the same /(C) as defined in Section [31 That is, if we temporarily forget the restriction 
to u; > 0, A^pcr equals — = — A'wcyi for uj G vanishes at the endpoints, and is 

periodic thereafter; this function jumps upward by 1 at each (positive) odd multiple of which 
we know to be the correct eigenvalues for these problems. (Another way of saying this sort of 
thing is that iVper is equal to the negative of N^^cyl i plus a series of unit step functions located 
at the eigenvalues.) At = this function equals 0, so the complete A'pcr is continuous there. 
For consistency, iVbdry must turn out to be when / + r is odd (and it does). 
If / + r is even. 



TT 

n=l 



so the jumps occur at the (positive) integer multiples of ^, as they should. The only complication 
is at w = 0, where this function approaches ^ from the right. Since the complete A^per is for 
UJ < 0, there is only half a step function at 0. In the full N, there should be a complete unit 
jump in the Neumann problem {I = 0) and no jump at all in the Dirichlet problem (l = 1). 
These corrections must come, of course, from the boundary term. 
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For it we have 

(-1)' 



bdry 



n=0 



2Luj{n+l) 



smu 



du 



[-Vf °° r2Luj{n+l) 



TT 



n=0 



sinu 



2Lu)n 



smu 



-1)' 



( /"OO 



2Lujn 



U 



TT 



sinn , (— 1)' 

du = -^r— 

u 2 



(see, e.g., [Ml 3.721.1]). This holds for lo > 0; thus A^'bdry heals the discontinuity at u; = 
in A'per when I = 1 and strengthens it to a unit jump when / = 0, as expected. 

Also noteworthy, although not unexpected, are that A^'bdry is independent of L - it is as- 
sociated with the physics of the boundary, not the finiteness of the region - and that it is 
nonoscillatory as a function of uj. Thus for N (and other global quantities) it seems proper to 
write 



Nav = AVcyl + A^bdry, 



-Aqsc — -Aperj 



in contrast to the definition of a^v and dose- 

If we kept only the contribution from the short orbits, we would get 



iVbd 



■1)' 



I r2Lu) 



TT 



smti 



u 



du, 



a fair approximation to the correct step function when Lu) is large. For comparison, our for- 
mula for /Obdry on the half-line also yields a step function at the origin of half that magnitude, 

H{uj), as one should expect [10] for only a single boundary. 



Now consider the regularized vacuum energy 
E{t)- 



d 1 
dt 2 



poo 

/ P('^)^~'^* du) = E'wcyl + EpcT + -E'bdry, 
JO 



where 



, , L d I 



L 



2iT dt t 27rt2 ' 



and 



Eperit) = ---n V(-l)"('+'^) / cos(2..nL)e--* 
dt Jo 



L d 
vr dt 



(l+r) 



n=l 



t2 + (2nL)2 



duj 

_l_d_ 

2tt di 



(l+r) 



t/2L 



n=l 



(t/2L)2 + ?i2- 



These series can be evaluated by [Ml 1.217.1,2] or their generalizations [Ml 1.445.2,3]: In the 
even case, 



Epcr {t) 



1 d 
27r dt 



TT 



coth 



7rt\ L 
2l) ~1 



TT 

8L 



of Trt\ 

csch — — 



L 



\2Lj 27rt2' 



In the odd case, 
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Expand these in power series in t, using [3i\ 1.411.8,12]: 



cschz = - - - + 0(z^ 
z 6 



cothz = - + - + 0(z^ 
z 3 



The terms of order t ^ cancel, and we get 



Epcr (t) 



241 

TT 

48L 



+ 0{t^) if / + r is even, 
+ 0{t^) if Z + r is odd. 



The first of these gives the well known renormalized vacuum energy (j2.7p for the one-dimensional 
Dirichlet problem [l = r = 1). It is also correct for the Neumann case, when any energy 
associated with the indiscretely quantized zero mode is neglected. The extra factor of — i in 



the mixed case is just like that for the antiperiodic case in Section [3] and [52] and similar to 
the factor of — | in the mixed case in three-dimensional electromagnetism [8]; this family of 
formulas can be obtained at the eigenfunction level by doubling the interval and removing the 
contribution of the even modes [18j . 

It remains to investigate the boundary energy. As usual it will be zero in the odd case. For 
the even case we have 



bdry 



2tt dt ^ 



n=0 



[sin(2Lo'L(n + 1)) — sin(2ajLn)] 



-LOt 



duj 



to 



27r ^ 



n=0 

oo 



27r ^ 



n=0 



[sin(2wL(n + 1)) - sin(2a;Ln)]e"'^* dco 
2L{n + 1) 2Ln 



This series is conditionally convergent and telescopes to zero. Alternatively, put the two terms 
in the summand over a common denominator: 



bdry 



(0 



~ 2Lf - {2Lfn{n + 1) 



The term for n = is 
(-1)' L 



i-iy 



In the other terms it is legitimate to expand the denominators before summing: 



(-1)' ^ 2Lt2 - (2L)^n(n + l) 



27r 



E 

n=l 



(2L)4n2(n+ 1) 



1 



(2Ln)' 



+ 



1 



(2L(n + l))2 



+ 



(-1/ 



AttL ^ n(n + 1) 

n=l 



4ttL 



+ 0(t'), 



the numerical sum being a textbook example [581 pp. 612-613] of a telescoping series that 
converges to 1. So all these terms exactly cancel the n = term in the limit of small t, and 
£^bdry(0) = 0. This result was not entirely obvious, since one might expect boundary energies 
(possibly infinite) at both ends with the same sign. 
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In the approximation of keeping only the short orbits, one gets a nonzero (and L-dependent) 
result, -EbdrylO) = ■ For the half-line we have 



2 dt Jo 4 ' ' 2 JO 

= 0. 

Now we turn to local quantities. First, integrate a to get a local analogue of the counting 
function. (This is the inverse Laplace transform of the diagonal value of the heat kernel; it is 
the quantity called /iOO in [25].) As expected. 



(Tav = -, j 0-per '^'^ = sin(2nLu;) = — iVper. 



n=l 

The boundary term is 



a,,,ycL; = - ^ (-1) . 

n=-oo ^ ' 

Integration of it over x yields the same A^bdry found before. The analogous calculations for the 
half-line give 



/ 



, (-1)' sin(2u;x) 

CJbdry AijJ 



27r X 

which is precisely the n = term in the sum above, and 



2-K 







in agreement with our previous result (14. ip for this case. 

At this point let's pause to compare the results with those from a direct summation of the 
eigenfunction expansion, specializing to the doubly Dirichlet case, I = r = \. The diagonal value 
of the (integrated) spectral kernel is 

inj/L)<uj j=l ^ ^ 

which reduces after some calculation to 

i=-[wL/7rJ 



The first two terms are precisely the step function -^(Aav + Apcr). The sum in the final term 
is the Dirichlet kernel introduced in any rigorous textbook on Fourier series (e.g., [38^ p. 22]). 
Thus 
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When 00 is large, this function develops sharp peaks near x = and x = L, in keeping with the 
boundary-effect picture we have had in mind all along. But instead of performing this sum, we 
can relate it to the closed-orbit calculation by the Poisson summation formula: 



2L a bdry 



luiL/n\ 

E » 

jr=-[wL/7rJ 



2-Kijx/L 



^2TTijx/L 



9{ujL 



vr 



l= — C< 

oo 

E 



[wL/ttJ 
[wL/tt] 



^2TTij{n+x/L) 



dj 



oo 

E 

i=—oo 



sin(2a;(x + nL)) 
X + nL 



Finally, we calculate the local energy density. The definition of the energy density of the 
massless scalar field in spatial dimension 1 in flat space is 



dx ^ 



dx , ~'~ 



dx"^ 



where ^ is called the conformal coupling parameter; different values of ^ correspond to different 
theories of the coupling of the field to gravity, but in flat space they are physically and mathe- 
matically equivalent apart from this one definition. The term multiplied by —4,^ equals 



d_ 

dx 



dx 



which classically vanishes upon integration by parts in free space or under either Dirichlet or 
pure Neumann boundary conditions. In dimension 1 there are only two distinguished values of ^, 
namely (the conformal, or Yamabe, choice) and j (which is needed for energy conservation in 
the presence of general boundary conditions if singular surface terms are to be avoided |46t E31 
126]). The energy density for general ^ is a convex combination of these two special cases. In the 
conformal case in dimension 1 the contribution of each normal mode to the vacuum expectation 
value is independent of x, so we know that Tqq is just E/L. (In particular, there is no boundary 
contribution.) In the case C = \ the contribution of the space derivatives is identical to that of 
the time derivatives, so we can write 



Too(i,x) = E{t,x) 



ld_ 

2 dt 



a{uj)e du; = £'wcyl(t) + £^per(i) + EhdTy{t,x), 



which now is indeed the energy formula with the integration stripped off. 
Clearly, for the first two terms we get the same old result. 



-^Wcyl(O) + £'pcr(0) 

The interesting term is 
Ehdry{t,x) = 



E 

T 



-ly d 



277 



oo 

di ^ 



oo „ 
(-1) 



2tt dt 

{-ly d I i\n.(i+r) 



cos(2l^(x + nLy)e~'^^ duj 



t2 + 4(x + nL)2 



-11 



27r 



oo 

E 

i=— oo 



(-1) 



n{l+r) 



t^ - 4(x + nLf 
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These sums do not appear in [3l] , but Mathematica |60j evaluates them in terms of hyperbolic 
functions of complex argument. Some hindsight reveals what is going on: The summations we 
used to evaluate the total energy (effectively the present sums with x = 0) can be written 

— cothlvra) = > rr, — cschlvra = > ^. 

a ^ -n? + a ^ 'n? + a? 

n=—oo 71=— oo 

Factoring the terms on the right, we see that these are just Mittag-Leffler expansions of the 
hyperbolic functions in terms of simple poles: 



27rcoth(7ra) = f — \ — ) 

^-^ \a + in a — in J 



and a similar formula for csch. So the summands with a general quadratic in the denominator 
can be treated by displacing the argument of the hyperbolic functions and letting n i— > —n in 
half the terms: 



oo ^ 

J2 (n + b)^ + a2 " + ^^)) + coth(7r(a - ib))], 



n=—oo 
oo 



^ f — 1)" TT 

J2 (n + 6)2 + a2 = ^[csch(7r(a + ib)) + csch(7r(a - ib))]. 

n=—oo ^ ' 

(These formulas are actually improvements on Mathematical s output, though equivalent.) 
Thus when / + r is even, 



^ , , (-1)' d 
£^bdry(t,x)--^- 



, / vrt i7rx\ , / nt ittx 

coth \ + coth 

\2L L J \2L L 



and when Z + r is odd, the formula is the same with coth replaced by csch. Differentiating first, 
and then working out the complex arithmetic, one gets in the even case 

E u ^^ _ (-1)'^ ^i^h^ (Ix) ^o^' (f ) - -o^^' (It) (f) 

[,i,l,2 + ^^^^2 ^.^2 (^)] 2 

and in the odd case 

-IYtt cosh (|1) cos (f ) [sinh^ (|1) - sin^ (^) + 2sinh2 (^) sin^ {^)] 



Ebdrj{t,x) 



1 2 



[sinh^ (|1) cos2 {^) + cosh^ (^) sin^ (f^)] ■ 
Specializing to / = 1, let us examine the leading terms at small t and small x. For the even 



case, 



-Ebdry(0,x) = ^ CSC — 

This is the renormalized boundary energy density in (j2.15p (Fig. [T]) . Its integral over the whole 
interval is not even finite, much less zero as formally expected. Near x = we have 

(and a corresponding expansion near x = L). On the other hand, if we expand in x first we get 

Ehdry{t,x) = - ^ CSch^ j +0(x^). 
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the remainder term being nonuniform in t. At a; = and small t this becomes 

£^bdry(i,0) ^ - 

- that is, infinitely negative! Plotting the exact £^bdry(^) ^) for various small but nonzero values of 
t reveals a steep rise as x approaches the boundary, followed by an even steeper plunge to negative 
values still closer to the boundary (cf. Fig. [2]). This behavior assures that E]^dry{t, x) dx = 
for any nonzero t, so that the total energy (j2.7p is independent of ^, as it must be. (The integral 
has been evaluated directly in the half- line case (14. 3p . and with the cutoff in place there is no 
obstacle to integrating the original spectral sum term by term and observing that the total- 
derivative terms integrate to zero.) This mathematical phenomenon was pointed out by Ford 
and Svaiter 



The situation for the odd case is very similar. Sbdry(i)a^) is now an odd function of 



X 



L 
2 ' 

with the same qualitative behavior as in the even case near x = and the inverted behavior at 
the other end of the interval. The formulas parallel to those above are 

TT / VrX \ / TTX 
rnt I 1 rsp I 



^^bdry(0,x) = ^ cot — j CSC — 



TT 



TTt\ , / vri \ 



-£^bdry(i, ^) = - ^ [jlj ^^'^^ [jlj + '^(^ )' 

-Ebdry(i,0) W - 

For the problem with a single boundary at x = we have 

Ehdry{t, ^) = ~ 2 at / ^bdryCu^, x)e"'^* '^^ ^ ~ 2 9t / TT cos(2wx)e"'^* duj 

{-ly d t _ {-ly t^-Ax'^ 

2tt di t2 + 4x2 - 27r (t2 + 4x2)2' 

Restricted to < x < L, this is exactly the n = term in the sum for the problem with two 
boundaries. (The other short path, n = —1, naturally gives a symmetrical contribution localized 
at the other boundary.) The leading terms are the same as found above for two boundaries. In 
particular, 

i-iy-^ 

-Bbdry(0,x) = g^^2 ' 

in agreement with the renormalized vacuum energy found by Romeo and Saharian [531 (3.21)]. 
At small X and fixed t we have 



-£'bdry(*> x) 



(-1) 

2tt 



1 12x2 ^/x^ 
+ 



t2 ' 



again giving some insight into the sharp spike of opposite sign that keeps the total boundary 
energy equal to zero as long as the regularization has not been removed. Indeed, in this case 
the integral of the exact function is elementary: 



t2 -4x2 X 

■ ax 



(t2 + 4x2)2 -^-^ t2 + 4^2 

as reported in Section [21 



= 0, (4.3) 

x=0 
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5 Conclusions 

One of our major concerns in working out these elementary models in complete detail has been 
to appraise the stationary-phase approximation (which here means completely discarding the 
nonperiodic closed orbits) and the approximation of ignoring all nonperiodic orbits except the 
shortest ones (which inevitably leads to a comparison with the problem on the half- line). Let 
us summarize the observations, in reverse order. 

For the local vacuum energy density, the contribution of a short orbit is equal to the corre- 
sponding term in the half- line problem (restricted to the interval, of course). These are good 
approximations to the "boundary" part of the exact answer near the corresponding boundary, 
which is the only place where they are large. These boundary terms are not zero or generally 
small, so the stationary-phase approximation is not very good here. 

For the local spectral density, or its unnamed indefinite integral, again the short orbit's 
contribution is equal to the half- line expression. Comparison with the Dirichlet kernel again 
suggests that these are good approximations near the boundary, and that these boundary terms 
are not very important elsewhere - although in this case they decay by virtue of oscillation 
(distributionally) more than by decrease of magnitude. 

For the total vacuum energy, globally renormalized, we found that the boundary terms in 
both the exact answer and the half-line expression are equal to zero; that is, the stationary-phase 
approximation is exact! The short orbits, however, give a nonzero result; it is of roughly the 
same order of magnitude as the exact Casimir energy coming from the periodic orbits, so it 
must be rejected as wrong. The apparent contradiction with the local result is explained by the 
difference between local and global renormalization. 

For the eigenvalue density, or its integral the counting function, the half-line theory and 
the exact theory agree (when both endpoints are counted), and they give the correction to the 
stationary-phase theory that is necessary to take account of the effect of the boundary conditions 
on the position of the lowest eigenvalue. The short-orbit expression differs, but it seems to be 
a reasonable approximation. 
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